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^  Pol arimetry  and  spectropolarimetry  are  optical  measurement  techniques  which  use  pol¬ 
arized  light  to  obtain  electrooptical  material  property  information.  These  tech¬ 
niques  are  being  used  to  make  measurements  on  infrared  electrooptical  materials. 
Infrared  materials  of  interest  are  those  which  find  use  in  two  dimensional  modula¬ 
tors,  i.e.  optical  processing  applications  and  infrared  scene  projectors.  Polarim-* 
etry  measurements  are  made  at  one  infrared  wavelength  at  a  time  using  laser  sources. 
Spectropol arimetry  measurements  are  made  over  an  entire  infrared  spectral  region, 
e.g.  2  to  14  urn.  A  Fourier  transform  infrared  spectrometer  is  the  source  of  radia¬ 
tion  in  this  case.  The  Mueller  matrix  formulation  can  be  used  to  describe  the  pol¬ 
arized  light  and  its  interaction  with  the  sample  and  measurement  system  optics.  A 
Mueller  matrix  of  the  sample  can  then  be  obtained  from  a  series  of  measurements  with 
different  input  polarization  states.  The  sample  Mueller  matrix  contains  information 
on  the  polarization  properties.  Electric  fields  are  imposed  on  the  sample  as  optical 
measurements  are  made.  Fundamental  constants  associated  with  the  sample  material 
can  be  derived  from  the  measured  Mueller  matrix  elements.  - 
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ABSTRACT 


Polarimetry  and  spectropolarimetry  are  optical  measurement  techniques  which  use  polarized  light  to  obtain 
electrooptical  material  property  information.  These  techniques  are  being  used  to  make  measurements  on 
infrared  electrooptical  materials.  Infrared  materials  of  interest  are  those  which  find  use  in  two  dimensional  mod¬ 
ulators,  i.e.  optical  processing  applications  and  infrared  scene  projectors.  Polarimetry  measurements  are  made  at 
one  infrared  wavelength  at  a  time  using  laser  sources.  Spectropolarimetry  measurements  are  made  over  an 
entire  infrared  spectral  region,  e.g.  2  to  14  pm.  A  Fourier  transform  infrared  spectrometer  is  the  source  of  radi¬ 
ation  in  this  case.  The  Mueller  matrix  formulation  can  be  used  to  describe  the  polarized  light  and  its  interaction 
with  the  sample  and  measurement  system  optics.  A  Mueller  matrix  of  the  sample  can  then  be  obtained  from  a 
series  of  measurements  with  different  input  polarization  states.  The  sample  Mueller  matrix  contains  information 
on  the  polarization  properties.  Electric  fields  are  imposed  on  the  sample  as  optical  measurements  are  made. 
Fundamental  constants  associated  with  the  sample  material  can  be  derived  from  the  measured  Mueller  matrix 
elements. 

1.  INTRODUCTION 

Bulk  properties  of  transmissive  materials  can  be  measured  using  the  experimental  methods  of  polarimetry  and 
spectropolarimetry.  We  show  in  this  paper  how  these  techniques  are  implemented  and  how  fundamental  proper¬ 
ties  may  be  derived  from  the  measurements. 

The  Mueller  matrix  formalism  is  used  to  process  the  polarimetric  information.  A  mathematical  background  on 
the  Mueller  formalism  and  crystal  physics  are  presented.  We  then  describe  details  of  an  infrared  laser  polarime- 
ter  and  an  infrared  spectropolarimeter.  Derivation  of  electrooptic  tensor  coefficients  directly  from  the  measured 
Mueller  matrix  is  discussed. 
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2.  MATHEMATICAL  BACKGROUND 

The  mathematics  for  describing  the  polarization  modulation  when  a  beam  of  optical  radiation  is  passed 
through  a  material  is  presented  here.  The  mathematical  formalism  for  polarized  light  is  discussed  first,  followed 
by  the  mathematics  necessary  to  describe  crystal  physics. 

Two  algebraic  systems  have  been  developed  for  the  solution  of  polarization  problems  in  optics,  the  Jones  for¬ 
malism  and  the  Mueller  formalism.  The  Jones  formalism  is  a  natural  consequence  of  the  mathematical  phase 
and  amplitude  description  of  light.  The  Mueller  formalism  comes  from  experimental  considerations  of  the 
intensity  measurements  of  polarized  light.  The  Mueller  formalism  is  used  for  the  present  work  and  is  described 
below.  Other  treatments  of  the  Mueller  and  Jones  formalisms  can  be  found  in  Gerrard  and  Burch^,  Azzam  and 
Bashara^,  Clarke  and  Grainger^,  and  Theocaris  and  Gdoutos'*. 
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2.1  The  Mueller  formalism 


The  Mueller  formalism  owes  its  name  to  Hans  Mueller,  who  built  upon  the  work  of  Stokes®,  Soleillet®,  and 
Perri  n  7  to  formalize  polarization  calculations  based  on  intensity.  This  work  was  also  done  during  the  1940’s  but 
originally  appeared  in  a  now  declassified  report  (Mueller®)  and  in  a  course  of  lectures  at  MIT  in  1945-1946. 

The  Mueller  formalism  uses  the  Stokes  vector,  a  four  element  real  vector,  to  represent  the  polarization  state. 
This  vector  is  often  represented  as 


S  = 


(1) 


where  s  0is  related  to  intensity,  s  i  to  preference  for  horizontal  polarization,  s  2to  preference  for  45°  polarization, 
and  s  3to  preference  for  right  circular  polarization.  The  inequality 


s0>s,  +  s,  +  s: 


(2) 


holds,  indicating  that  some  part  of  the  light  intensity  may  be  in  a  depolarized  state. 

The  Mueller  matrix  is  a  4x4  matrix  representing  optical  elements  or  an  optical  system,  i.e. 

f  mu  m-12  m13 

m  21  mzz  m  23 

M  - 

m3l  m32  m33 

\m41  m42  m43 


(3) 


where  the  m  <,  are  real  numbers. 

The  Stokes  vector  describes  completely  the  polarization  state  of  the  light  in  terms  of  intensity.  The  Mueller 
matrix  is  an  empirical  entity  for  a  given  ray  path  through  an  optical  element  or  system,  taking  the  form  required 
to  produce  the  correct  output  polarization  state.  There  is  redundancy  built  into  the  Mueller  matrix,  since  only 
seven  of  its  elements  are  independent^  if  there  is  no  depolarization  in  the  optical  system.  In  the  most  general 
case,  the  Mueller  matrix  can  have  sixteen  independent  elements  1®,  however  not  every  4x4  Mueller  matrix  is  a 
physically  realizable  polarizing  element.  The  latter  aspect  of  the  Mueller  formalism  is  discussed  by  Cloude^^»^^. 

Matrix  products  represent  the  action  of  sequential  optical  elements  on  a  light  beam.  In  the  Mueller  formalism 
this  is  represented  by 

i 

A'f  s  =  A  /  K. . .  .W  ZM  ,  =  ]””[  M  k  ^4) 

k-K 

where  M  5  is  the  system  Mueller  matrix  and  the  M  k  are  the  element  matrices.  The  element  matrices  are 

multiplied  on  the  left  in  the  order  that  the  light  encounters  the  elements.  The  Mueller  matrices  are  in  general 
functions  of  wavelength  and  their  own  orientation  ©with  respect  to  some  arbitrary  reference,  i.e.  M  =  M(\ ,  0). 


The  Mueller  matrix  can,  unlike  the  Jones  matrix,  represent  a  depolarizer  or  scatterer,  and  van  de  Hulst^ 
describes  this  property.  The  Mueller  matrix  formalism  has  two  advantages  for  experimental  work  over  the  Jones 
formalism.  The  intensity  is  represented  explicitly  in  the  Mueller  formalism,  and  scattering  can  be  included  in  the 
calculations. 

2.2  Crystal  physics 

In  general,  crystals  are  anisotropic  with  respect  to  the  propagation  of  electromagnetic  waves  in  the  presence  of 
static  electric  fields.  The  relationship  between  the  displacement  and  the  field  is  (MKS  units) 

D,  =  z„E,  (5) 


where  e  ,t  is  the  dielectric  or  permittivity  tensor^.  The  impermeability  tensor  is  defined  as 


n,;  =  e0(e"  )„ 


(6) 


where  e " 1  is  the  inverse  of  the  dielectric  tensor  and  e  0  is  the  permittivity  of  free  space.  The  principal  indices  of 
refraction  n„ny,  and  n  .  are  related  to  the  principal  values  of  the  impermeability  tensor  and  the  principal  values 
of  the  permittivity  tensor  by 

1  e o  (?) 

,2  n‘*  f  ' 

l  -  kk 

The  electrooptic  coefficients  are  defined  by  the  equation 

nI/(F)-niy(0)-Aii1/  =  rl/t£t-si/U£,£(  +  O(£'’)  (8) 


where  q,;  is  a  function  of  the  applied  field  E,  rl/k  are  the  linear  electrooptic  tensor  or  Pockels  electrooptic 
coefficients,  and  the  s  w*,are  the  quadratic  or  Kerr  electrooptic  coefficients.  Terms  higher  than  quadratic  are 
small  and  are  neglected. 

The  linear  electrooptic  tensor  is  of  third  rank  with  3^  elements;  however,  symmetiy  enables  a  contraction.  If 
the  medium  is  lossless  and  optically  inactive,  e is  a  symmetric  tensor,  n„  >s  a  symmetric  tensor,  rl/k  has  symme¬ 
tiy  where  coefficients  with  permuted  first  and  second  indices  are  equal,  and  s  ,,*,has  symmetry  where  coefficients 
with  permuted  first  and  second  indices  are  equal  and  coefficients  with  third  and  fourth  coefficients  are  equal. 

The  linear  electrooptic  coefficients  are  assigned  two  indices  so  that  they  are  r  where  l  runs  from  1  to  6  and  k 

runs  from  1  to  3.  The  quadratic  coefficients  are  assigned  two  indices  so  that  they  become  s ,,  where  :  runs  from  1 

to  6  and,/  runs  from  1  to  6.  For  a  given  crystal  symmetiy  class,  the  form  of  the  electrooptic  tenser  is  available 
from  standard  tables^. 

Light  propagating  in  anisotropic  materials  experiences  refractive  index  and  phase  velocity  dependent  upon  the 
propagation  direction  and  polarization  state.  The  refractive  index  for  propagation  in  some  direction 
k  =  xi  +  y  j  +  zk  can  be  obtained  from  the  index  ellipsoid.  In  the  principal  coordinate  system  the  index  ellipsoid 
is  given  by 


(9) 


in  the  absence  of  an  applied  electric  field.  The  length  of  the  semimajor  and  semiminor  axes  of  the  ellipse  formed 
by  the  intersection  of  this  index  ellipsoid  and  a  plane  normal  to  the  propagation  direction  (both  geometric  objects 
being  centered  at  the  origin)  are  the  two  principal  indices  of  refraction  for  that  propagation  direction.  If  two  of 
the  three  indices  of  the  index  ellipsoid  are  equal,  the  crystal  is  defined  to  be  uniaxial  and  the  equation  for  the 
index  ellipsoid  is 


-V 


(10) 


Where  there  are  three  distinct  principal  indices,  the  crystal  is  defined  as  biaxial. 

The  alteration  of  the  index  ellipsoid  of  a  crystal  with  application  of  an  electric  and/or  magnetic  field  is  useful 
for  the  modulation  of  polarization  state.  The  equation  for  the  index  ellipsoid  of  a  crystal  in  an  electric  field  is 

ri(/(f  1  (H) 


or 


(n„(0)+  At1|/).v,.v/=  1 


This  equation  can  be  written  as 


1  .  (  1 


njt 
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a  -  +y‘\-a  +  *\z 


l  .  ( i 


•y 
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+  2yz  A  -  +2.vz  A  -  +  2.vy  A  -  =1 
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+  z2 1  — r  +  A(  — 


or 


ru^*  +  si*^k  +  2sMF2F  3*2st5F3F  |  +  2s  l6F  j  F  2 


y2  2  +  r2*£'*  +  s2tFt  +  2s245'2F3-'-2s2SF3F  1  +  2s26F  ,F; 


+  z  ^2  +  r3t^»  +  s3*^*  +  2s34F2F3-*-2S;,5F3F  |+2s36F|F; 

+  2yz(r4tFts4nF4  +  2s44F2F3  +  2s4SF3F,  *2 s46F,F2) 

+  2zx(r5tFt  +  s5(tF2  +  2sssF 3F  ,  +  2 s56F,F2) 

+  2xy  (r64  Ft  +  s  bk  E  k  +  2s64F2F3*  2  s65  F  3F  ,  +  2  s  66  F  ,F  2)  =  1 


(12) 


(13) 


(14) 


where  the  F  *  are  components  of  the  electric  field  along  the  principal  axes  and  repeated  indices  are  summed. 


If  the  quadratic  coefficients  are  assumed  to  be  small  and  only  the  linear  coefficients  are  retained,  then 


2  3 


A  n  Y.FltE  ‘ 


(15) 


and  k  =  1,2,3  corresponds  to  the  principal  axes  x,  y,  and  z.  For  a  uniaxial  crystal  in  a  field  along  z,  the  index 
ellipsoid  becomes 


X 2  V 2  2  2 

—  +  V-  +  2r63£*xy=  I- 

n,  na  n. 


(16) 


A  coordinate  rotation  can  be  done  to  obtain  the  major  axes  of  the  new  ellipsoid.  In  the  present  example,  this 
yields  the  new  ellipsoid 


-Vr63£\V2+[  -2-r63Ez)y'2  +  ^=  1 

ni  J  \n  i  J  n\ 


(17) 


The  new  and  old  z  axes  are  the  same,  but  the  x’  and  y’  axes  have  rotated  around  the  z  axis. 
The  refractive  indices  along  the  new  x  and  y  axes  are: 

1  3  r 

n*  =  n0--n0r63E  z 

1  3  , 

ny  =  n0+  ~n0rb3Ez. 


(18) 


F or  light  linearly  polarized  at  45  degrees,  the  x  and  y  components  experience  the  different  refractive  indices  n  x ' 
and  ny‘.  The  birefringence  is  defined  as  the  index  difference  ny'  -  nx'.  Since  the  phase  velocities  of  the  x  andy 
are  different,  there  is  a  phase  retardation  f  (in  radians)  between  the  x  and  y  components  of  E  given  by 


r  to  ,  ...  2n  3 

r  =  -(ny  -  nx  )/  =  —n0r63E zl 

C  A. 


(19) 


where  Z  is  the  path  length  of  light  in  the  crystal.  The  electric  field  of  the  incident  light  beam  is: 


(20) 


After  transmission  through  the  crystal,  the  electric  field  is: 


*  r  /  i  r  /  2  <>  <  - 1  r  /  2  \ 

=  t  ( G  X  +  Q  y  ) 


(21) 


v  2 


If  the  path  length  and  birefringence  are  selected  such  that  1'  =  n,  the  modulated  crystal  acts  as  a  half  wave 
linear  retarder  and  the  transmitted  light  has  field  components 


(22) 


-L  £ ( e ,n/2 .-C '  -  g "’,/2 y ' )  =  £ ( e'^2 .v  '  -  e y ' ) 

V  2  V  2 


The  axis  of  linear  polarization  of  the  incident  beam  has  been  rotated  by  90  degrees  by  the  phase  retardation  of  n 
radians  or  one-half  wavelength.  The  incident  linear  polarization  state  has  been  rotated  into  the  orthogonal 
polarization  state.  An  analyzer  at  the  output  end  of  the  crystal  aligned  with  the  incident  (or  unmodulated)  plane 
of  polarization  will  block  the  modulated  beam.  For  an  arbitrary  applied  voltage  producing  a  phase  retardation  of 
T  the  analyzer  transmits  a  fractional  intensity  cos2r.  This  is  the  principle  of  the  Pockels  cell. 

3.  POLARIMETRY 

Polarimetiy,  or  transmission  ellipsometry,  is  the  branch  of  ellipsometry  concerned  with  measuring  the  polar¬ 
ization  state  of  a  light  beam  and  deducing  the  polarizing  properties  of  bulk  materials.  A  polarimeter  is  an  optical 
instrument  used  for  the  determination  of  the  polarization  state  of  a  light  beam.  Given  the  knowledge  of  how  the 
polarimeter  itself  acts  on  the  light,  the  polarization  change  produced  by  inserting  a  bulk  sample  into  the  polari¬ 
meter  may  be  determined.  Implementation  of  polarimetry  can  be  identical  to  ellipsometry  except  that  the  optical 
path  goes  completely  through  the  sample. 

The  infrared  laser  polarimeter  described  here  is  designed  to  be  operated  over  the  3  to  14  pm  spectral  region 
using  laser  sources  at  several  wavelengths.  This  spectral  region  is  of  great  interest  for  evaluation  of  materials 
used  as  elements  of  optical  processing  systems  or  thermal  imaging  systems. 

The  optical  configuration  used  is  called  the  dual  rotating  compensator  ellipsometer  when  it  is  used  in  ellip¬ 
sometry.  The  data  processing  prescription  which  accompanies  these  optics  has  been  described  by  Azzam14  and 
by  Hauge1^  This  polarimeter  uses  two  fixed  polarizers  and  two  rotating  waveplates  in  a  configuration  discussed 
below.  Mueller  matrix  elements  and  Stokes  vectors  are  used  here  to  represent  the  polarization  elements  and 
polarized  light,  respectively.  The  Mueller  matrix  formulation  is  used  because  it  is  preferable  for  experimental 
work  where  scattering  and  depolarization  measurements  are  required.  Two  waveplates  are  rotated  at  different 
but  harmonic  rates  and  a  modulation  of  the  detected  intensity  results.  The  Mueller  matrix  of  the  sample  is  found 
through  a  relationship  between  the  Fourier  coefficients  of  a  series  representing  the  modulation  and  the  elements 
of  the  sample  matrix. 

Figure  1  shows  a  block  diagram  of  the  polarimeter.  The  polarimeter  has  five  sections,  the  laser  source,  the 
polarizing  optics,  the  sample,  the  analyzing  optics,  and  the  detector. 
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Figure  1:  Polarimeter  block  diagram 


The  polarizing  optics  consist  of  a  polarizer  and  a  quarter  wave  plate  both  mounted  in  computer  controlled 
rotating  stages.  The  sample  region  is  followed  by  the  analyzing  optics  which  consist  of  a  quarter  wave  plate  in  a 
computer  controlled  rotating  stage  followed  by  a  polarizer.  Polarization  sensitivity  of  the  detector  is  not  impor¬ 
tant  since  the  orientation  of  the  final  polarizer  does  not  change. 

Figure  2  shows  the  instrument  in  more  detail.  This  initial  work  has  been  done  with  a  CO2  laser.  The  beam  is 
chopped  at  about  1  kHz.  A  beamsplitter  is  used  to  divert  a  fraction  of  the  beam  to  monitor  the  power  output  of 
the  laser.  The  remaining  beam  then  passes  through  the  polarimeter  elements  including  the  sample,  and  is 
focused  with  a  ZnSe  lens  onto  a  HgCdTe  detector. 
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Figure  2:  Polarimeter  optics  diagram 


This  polarimeter  measures  a  chopped  signal  which  is  modulated  by  rotating  the  polarizing  optical  elements. 
The  elements  of  the  Mueller  matrix  can  be  encoded  on  the  modulated  signal.  The  output  signal  is  then  Fourier 
analyzed  to  determine  the  Mueller  matrix  elements.  This  particular  implementation,  due  to  Azzam^,  uses  two 
aligned  and  fixed  linear  polarizers  and  two  rotating  quarter  wave  retarders  as  shown  in  Figure  3. 
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Figure  3:  Polarizing  elements  and  rotation  rates 


The  second  retarder  is  rotated  at  a  rate  of  five  times  that  of  the  first.  This  generates  twelve  harmonic  frequen¬ 
cies  in  the  Fourier  spectrum  as  shown  below. 

The  Mueller  matrix  for  the  system  is 

L2  R2(0)  M  Ri(0)  Li  (23) 

where  L  indicates  a  linear  polarizer,  R  indicates  a  linear  retarder,  and  M  is  the  sample.  Mueller  matrices  are 
then  substituted  for:  a  linear  retarder  with  quarter  wave  retardation  and  a  fast  axis  at  0°  for  Rj  and  R2;  a 
horizontal  linear  polarizer  Lq;  a  linear  horizontal  polarizer  for  Lj;  and,  a  sample  for  M.  The  detected  intensity  is 
given  by 


1  =  cAMP 


(24) 


where  P  =  R!  LiS  is  the  Stokes  vector  of  light  leaving  the  polarizing  optics  with  S  being  the  Stokes  vector  of  the 
light  from  the  source,  A  =  L2  R2  is  the  Mueller  matrix  of  the  analyzing  optics,  M  is  the  Mueller  matrix  of  the 
sample,  and  c  is  a  proportionality  constant  obtained  from  the  absolute  intensity.  Explicitly, 
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/  -c  Y.  atP,mi, 
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(25) 


or 


4 

I  =  C 

i.j-i 


(26) 


where  the  a ,  are  the  elements  of  A,  the  p  ,  are  the  elements  of  P,  the  mly  are  the  elements  of  the  Mueller  matrix 
M,  and  where 


V-u  ~  aiPj- 


(27) 


The  order  of  the  multiplication  can  be  changed  as  above  because  we  are  only  measuring  the  intensity,  i.e.  the 
first  element  of  the  Stokes  vector.  Only  the  first  row  of  the  matrix  A  is  involved  in  the  calculation,  i.e. 
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The  p are  then  given  by 


Ull  - 


1 , 

412  =  cos220, 

4.3  =  sin29cos20, 

4.4  =  sin20, 

42,  =  cos2  109, 

4  22  =  cos2  20  cos  2  100, 

4 23  =  sin  20 cos 20 cos 2  109, 

4 24  =  sin  20cos2  100 , 

43,  =  sin  lO0cos  100, 

4  32  =  cos2  20  sin  100  cos  100, 

4  33  =  sin  20  cos  20  sin  100  cos  109, 

434  =  sin20sin  lO0coslOO, 

4  4,  =  -sin  100, 

4  42  =  -  cos2  20  sin  100, 

4  43  =  -  sin  20cos20sin  1  00 , 

4 44  =  ~  sin  20 sin  1 00,  (30) 

when  the  rotation  ratio  is  5:1.  These  equations  can  be  expanded  in  a  Fourier  series  to  yield  the  Fourier 
coefficients  which  are  functions  of  the  Mueller  matrix  elements. 

The  inversion  of  these  relations  gives  the  Mueller  matrix  elements  in  terms  of  the  Fourier  coefficients: 

m  |,  =ao-a2  +  a8-a,0  +  a,2, 

m,2  =  2a2-2a8-2a,2, 

n’t  1 3  =  2b2  +  2b„  -  2b  ,2 , 

m  ,  4  =  b  ,  -  2  b  , ,  =  b,  +  2b9  =  b,  +  b9-b,|, 

m2 ,  =  -2a8  +  2«,0-2a,2, 

m22  =  4 a 8  +  4 a  , 2 . 

m23  =  -4 b8  +  4b  ,2, 

m24  =  -4b 9  =  4b  ,,  =  2 (-b9  +  b  , , ). 

m3,  =  -  2b8  +  2b  ,0  -  2b  ,2  , 

m32  =  4b8  +  4b, 2, 

m33  =  4a8  -  4a 

m34  =  4a,  =  -4a  ,,  =  2(a„  -  a  , , ), 


(31) 


m4i  ~  2b3  b  s  —  b5  +  2b7-(b3  bs+b7), 

m42  =  -4b3  =  -4b7  =  -2(b3  +  b7), 
m43  =  -4a3  =  4a7  =  2(-a3  +  a7), 
m44  =  -2a4  =  2a6  =  (a6-a4). 

The  5:1  rotation  ratio  is  not  the  only  ratio  which  can  be  used  to  determine  Mueller  matrix  elements,  but  it  is  the 
lowest  ratio  where  the  expressions  for  the  Fourier  coefficients  may  be  inverted. 

4.  SPECTROPOLARIMETRY 

Spectropolarimetry  is  the  spectroscopic  study  of  the  polarization  properties  of  materials.  Thi3  measurement 
technique  is  a  generalization  of  conventional  optical  spectroscopy  and  of  conventional  polarimetry.  Commercial 
spectropolarimeters  and  spectropolarimeters  in  the  literature  are  usually  concerned  with  instruments  which  pass 
one  wavelength  at  a  time  through  a  sample.  These  instruments  typically  use  a  scanning  monochrometer  source. 
Because  of  the  polarization  effects  that  arise  from  diffraction  gratings,  and  the  length  of  time  that  would  be  noces- 
sary  to  obtain  sample  Mueller  matrices  at  each  wavelength,  these  spectropolarimeters  are  unsatisfactory  for 
acquisition  of  precision  polarimetric  data.  A  spectropolarimeter  based  on  a  Fourier  transform  infrared  (FTIR) 
spectrometer  has  been  designed  and  built!6,l /.  This  instrument  retains  the  advantages  of  a  FTIR  spectrometer 
over  conventional  spectrometers  and  is  uniquely  suited  to  collect  polarimetric  data  quickly  over  large  wavelength 
ranges.  An  optical  diagram  of  tne  FTIR  spectrometer  as  modified  to  be  a  spectropolarimeter  is  shown  in  Figure 
4.  This  spectropolarimeter  contains  the  same  polarization  element  configuration  used  in  the  laser  polarimeter 
described  above.  Spectrometer  data  at  each  orientation  of  the  polarization  elements  as  described  for  the  polari¬ 
meter  is  collected  and  integrated  to  yield  the  Mueller  matrix  of  a  sample  as  a  function  of  wavelength. 


Figure  4:  Spectropclarimeter  optics 


/ 

5.  MEASUREMENT  OF  ELECTROOPTIC  COEFFICIENTS 


The  goal  of  this  work  is  to  obtain  electrooptic  coefficients  from  Mueller  matrices.  The  method  by  which  this  is 
accomplished  in  described  in  this  section. 

The  application  of  an  electric  field  across  a  crystal  produces  an  index  change  as  described.  Specific  index 
changes  for  a  sample  material  is  best  presented  in  tabular  form.  Table  1  gives  the  index  ellipsoid  and  the  indices 
along  the  new  principal  axes.  The  new  principal  indices  are  obtained  by  solving  an  eigenvalue  problem.  For 
example,  for  a  cubic  material  with  a  field  perpendicular  to  the  (111 )  plane,  the  index  ellipsoid  is 
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and  the  eigenvalue  problem  is 


2r.,£  2/ 


The  secular  equation  is  then 
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and  the  roots  of  this  equation  are  the  new  principal  indices. 

The  phase  retardation  accumulated  by  polarized  light  in  traversing  a  medium  with  anisotropic  properties  is 

given  by 


r  =  2it(an-ni))Z./\ 


(35) 


where  L  is  the  medium  thickness  in  the  direction  of  propagation,  \  is  the  wavelength  of  light,  and  na  ,  nbare  the 
indices  experienced  in  two  orthogonal  directions  perpendicular  to  the  direction  of  propagation.  In  the  longitudi¬ 
nal  mode  of  operation,  the  electric  field  and  propagation  direction  are  both  along  the  z  axis.  The  refractive  indices 
experienced  by  the  light  are  in  the  plane  containing  the  x  and  y  principal  axes.  If  the  light  polarization  and 
crystal  are  aligned  so  that  the  polarization  is  45°  from  either  principal  axis,  the  phase  retardation  will  be 

r  =  2rt(n./-n/)/./\  (36) 


where  n.,' ,  nx‘  are  the  (new)  principal  indices  with  the  field  applied.  (For  crystals  with  natural  birefringence  and 
with  no  electric  field,  these  indices  may  just  be  the  principal  indices.) 


The  phase  delays  for  light  polarized  at  45°  to  the  principal  axes  of  the  material  of  Table  1  can  now  be  calcu¬ 
lated.  The  principal  indices  of  the  materials  for  an  electric  field  parallel  to  the  £and  optical  axes  (longitudinal 
mode)  are  given  in  the  first  section  of  the  table.  The  phase  retardation  for  the  43m  cubic  material  is 

F  cubic  =  2nn„r  41  E  L/k.  <37) 


If  the  electric  field  is  expressed  in  terms  of  electric  potential  and  charge  separation,  i.e.  E  =  V/d,  then  the  phase 
retardation  is 


r:Z<  =  2itn03r41l//\  (38) 

since  the  charge  separation  d  is  equal  to  the  optical  path  through  the  crystal  L.  The  phase  retardation  is  linear 
with  applied  voltage.  A  retardation  of  Jt  in  the  longitudinal  mode  will  require  a  voltage 
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(39) 


The  half  wave  voltage  for  CdTe  is  about  44  kV  and  for  GaAs  it  is  about  98  kV  for  10.6  pm  radiation.  Since  the 
voltage  is  directly  proportional  to  wavelength,  the  half  wave  voltage  would  be  considerably  smaller  at  visible 
wavelengths. 

The  principal  indices  for  the  transverse  mode  are  given  in  the  second  section  of  the  table.  The  phase  retarda¬ 
tion  for  4  3  m  cubic  material  is 

rtus,(  =  2n  nlr^EL/X.  (40) 


This  is  the  same  expression  as  for  the  longitudinal  mode.  In  this  mode,  the  charge  separation  is  not  the  same  as 
the  optical  path,  so  that  when  E  is  given  as  V/d,  the  phase  delay  is  given  as 

r  ‘cubic  =  2un3r  AlV  L/dK  <41> 


The  half  wave  voltage  is  given  by 
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(42) 


The  half  wave  voltages  for  CdTe  and  GaAs  for  an  optical  path  of  50mm  and  electrode  separation  of  5mm  are 
4.4kV  and  9.8kV  respectively. 

The  Mueller  matrix  formalism  representation  of  a  retarder  with  fast  axis  at  arbitrary  orientation  angle  0 is 


/ 1  0  0 
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(43) 


where  the  retardance  is  &  If  the  retarder  fast  axis  is  assumed  to  be  at  0°,  the  matrix  becomes,  substituting  for  6 
the  retardance  of  the  crystal, 


1  0  0 

0  1  0 

0  0  cos  —  n3 r  (]  l'  — 

K  cl 

2  jt  L 

0  0  -  sin  —  n3i\X[  I'  — 

A.  Cl 

It  is  now  clear  that  the  electrooptic  coefficient  r  4 ,  can  be  obtained  from  the  measured  Mueller  matrix. 

Note  that  for  purposes  of  obtaining  the  electrooptic  coefficient  experimentally,  the  fast  axis  of  a  electrooptic 
crystal  acting  as  an  ideal  retarder  can  be  at  any  orientation.  The  (4,4)  matrix  element  of  the  matrix  for  a  retarder 
with  fast  axis  at  angle  0  is  independent  of  fast  axis  orientation,  and  the  fast  axis  orientation  can  be  eliminated 
elsewhere  by  adding  the  (2,2)  and  (3,3)  matrix  elements  or  squaring  and  adding  elements  in  the  fourth  row  and 
column.  Given  a  measured  Mueller  matrix  of  a  crystal,  a  known  applied  voltage,  and  a  known  refractive  index, 
the  electrooptic  coefficient  r  41  can  easily  be  obtained. 

6.  SUMMARY  AND  CONCLUSIONS 

Two  instruments,  an  infrared  laser  polarimeter  and  an  infrared  spectropolarimeter  have  been  described. 

Using  the  Mueller  formalism  for  the  treatment  of  polarization,  it  has  been  shown  that  fundamental  constants  of 
sample  electrooptic  materials  may  be  derived  from  measurement  procedures  described  in  this  paper.  Experimen¬ 
tal  procedures  must  be  accomplished  with  careful  preparation  because  of  the  many  difficulties  associated  with 
these  sensitive  polarization  measurements.  The  difficulty  of  obtaining  data  is  commensurate  with  the  difficulty 
of  obtaining  crystalline  samples,  and  the  equipment  investment  is  orders  of  magnitude  larger  than  the  already 
high  cost  of  the  crystalline  samples. 
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Table  1:  Electrooptical  properties  of  cubic  43mcrystals 
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